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Two ubiquitous features of frustrated spin systems stand out: massive degeneracy of their ground 
states and flat, or dispersionless, excitation branches. In real materials, the former is frequently 
lifted by secondary interactions or quantum fluctuations, in favor of an ordered or spin-liquid state, 
but the latter often survive. We demonstrate that flat modes may precipitate remarkably strong 
quantum effects even in the systems that are otherwise written off as almost entirely classical. The 
resultant spectral features should be reminiscent of the quasiparticle breakdown in quantum systems, 
only here the effect is strongly amplified by the flatness of spin-excitation branches, leading to the 
damping that is not vanishingly small even at S^l. We provide a theoretical analysis of excitation 
spectrum of the S = 5/2 iron-jarosite to illustrate our findings and to suggest further studies of this 
and other frustrated spin systems. 

PACS numbers: 75.10.Jm, 75.30.Ds, 75.50.Ee, 75.40.Gb, 78.70.Nx 


I. INTRODUCTION 

Ever since their inception in the 1950s,frustrated 
spin systems have been a source of new ideas for a wide 
variety of problems; unconventional superconductivity,^^ 
order-by-disorder phenomena, SI and correlated spin- 
liquid states, SI are among them. In the core of this fer¬ 
tility is the near degeneracy between a vast number of 
spin configurations, originating from competing interac¬ 
tions that are favoring mutually exclusive ground states. 
In frustrated magnetic materials and their models this 
massive degeneracy is responsible for an extreme sen¬ 
sitivity to subtle symmetry breaking effects,SEI strongly 
amplified role of subleading coupling terms, SI hierarchy of 
emergent energy scales, SI and order-by-disorder effects by 
thermal,Sni and quantum fluctuations.SSESl 

Concomitant of the ground-state degeneracy is an¬ 
other hallmark feature of the frustrat ed spin sys tems: 
flat excitation branches at low energiesThey 
owe their origin to both the topological structure of 
the underlying lattices that facilitate frustration and 
the insufficient constraint on the manifold of spin con¬ 
figurations. A subclass of frustrated magnets that 
exhibits flat mo des pr ominently is the kagome-lattice 

antiferromagnets.^^HIlMMI Uj^der the influence of sub¬ 
leading interactions, majority of the known kagome- 
lattice antiferromagnets order magnet ically with spins 
forming non-collinear structures^IMIHSTl often 

reminiscent of the classical 120° motif on each triangle, 
Fig. a). Such a pattern is also emblematic of the geo¬ 
metric frustration, manifesting a compromise reached by 
spins locally to partially satisfy their antiferromagnetic 
trends. 

The following aspect of this picture is crucial. The non- 
collinearity of the ordered spin pattern implies strong 
nonlinear, anharmonic effects.!^ The role of such effects 
in the ground-state selection of frustrated s ystems has 
been recognized since the early days of the fiel d^ l ^^ l ^ ^and, 
recently, an accurate, systematic treatments of the quan¬ 
tum order-by-disorder effect due to them has received 


significant attention. h^ l ^^ l 

On the other hand, their role in the excitation spectra 
of the kagome-lattice antiferromagnets has been hardly 
touched upon. In this work, we demonstrate that the 
nonlinear terms can be particularly important in the 
spectral properties of the flat-band frustrated magnets, 
leading to spectacularly strong quantum effects even in 
the systems that are assumed almost classical. The resul¬ 
tant spectral features bear a remarkable similarity to the 
quasiparticle breakdown signatures i n quan tum spin- and 
Bose-liquids, such as superfluid ^He,^332] exhibit 

characteristic termination points and ranges of energies 
where single-particle excitations are not well-defined and 
are dominated instead by broad continua. 

It is usually assumed that such drastic effects can only 
occur in the systems that are inherently quantum in 
nature.^ ^ l ^^ l ^^ l In our case, their origin is in the near reso¬ 
nance decay of the “normal”, i.e. dispersive, modes into 
pairs of the flat-mode excitations facilitated by the non¬ 
linear couplings. As such, the effect is strongly amplified 
by the density of states of the flat modes and is very sig¬ 
nificant even for large-spin systems that can otherwise 
appear as purely classical, resulting in the damping ef¬ 
fect Tk/^k^l- While in the following we give a detailed 
account of the spectral properties of a specific kagome- 
lattice antiferromagnet, S = 5/2 Fe-jarosite, encouraging 
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FIG. I: (Color online) (a) q = 0 type of spin ordering on the 
kagome lattice, (b) Directions of the DM vectors. Arrows on 
the bonds show the ordering of Si and Sj in 
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its further investigation by inelastic neutron scattering, 
the outlined scenario should be applicable to a wide v a- 
riety of other flat-band frustrated spin systems 

The paper is organized as follows. In Section [TT] we 
provide a general argument for the resonant-like decay 
to exist in the frustrated flat-band systems and lay out a 
qualitative expectation, which highlights an unusual phe¬ 
nomenon: decays remain significant even in the S ^ oc 
limit. In Section [III] we provide our results for the decay- 
dominated spectral features in Fe-jarosite. Section |IV| 
gives a brief summary. Technical details are delegated to 
Appendix 

II. NONLINEAR COUPLING AND 

RESONANT-LIKE DECAYS 

Because of the noncollinear structure of the ground- 
state spin configuration, the interacting spin excitations 
in the kagome-lattice antiferromagnets are described by 

H = + ^ Y. ^qp^k^tq&k^Mk + h-c.,(l) 

p+q=k 

where the first term accounts for the spin-wave energies 
while the second is an outcome of the anharmonic cou¬ 
pling of spins that results in the mutual transitions be¬ 
tween excitation branches, see Appendix [^ for details. 
Specifically, it couples dispersive excitations with the flat 
modes, allowing for the resonance-like decay of the former 
into the pairs of the latter. We note that this general form 
of bosonic Hamiltonian § occurs in a variety of con¬ 
texts, including o ther frustrated antiferromagnets with 
noncollinear order,l^^^^ as well as spin-liquids,^^ valence- 
bond solids ,1^ and Bose-liquids.^^ 

The full extend of the the 1 /^'-expansion also involves 
quartic and source cubic terms, see Refs. 1121281 Then, 
the magnon Green’s function for the branch /i is 

^/ik(^) — ^ ~ ^Aik “ 5 (2) 

in which the self-energy includes all such terms. 

However, it is only decay terms in 0 that are responsi¬ 
ble for the resonance-like decay phenomenon discussed in 
this work. Given the off-resonance character of the source 
term, the Hartree-Fock nature of the quart ic terms, and 
the large-5' limit of the problem, one can safely approx¬ 
imate the self-energy by its on-shell imaginary part, i.e. 
by I]^k(w)U=e^k ~ -*r^k- The decay rate T^k in the 
lowest-order approximation is given by 

r^k = I E I^Xq;kl'<^ - e.q - e,k-q), (3) 

q,zy?7 

where sum is over the branches of the decay products 
and an explicit form of the vertex ^q^k_q.k is given in 
Appendix]^ With that, evaluation of the spectral func¬ 
tion A^\^{uj) = —{l/7r)lmGf^\i{uj) is also straightforward. 


The anharmonic cubic terms appear in the Hamiltoni¬ 
ans of the noncollinear magnets via bosonization of the 
terms that have a form ^ 5f 5f^^^ in the local reference 
frame of the ordered moments.^ Because of that, cu¬ 
bic vertices in ^ necessarily scale with the spin value 
as ^q^p-k cx: v^5. Since the energies of the decay prod¬ 
ucts scale as oc 5, it follows that F^k in @ must be 
spin-independent. Contributions to the decay rate from 
the higher-order terms should then follow a natural 1/5 
expansion with the exception of some special contours 
in k-space where a log(5) enhancement in 0 is pro¬ 
duced due to va n Hove singularities of the two-magnon 
continuumTherefore, one can conclude that for 
magnets with 5^1, damping of higher-energy magnetic 
excitations due to decays into lower-energy ones must be 
small compared to the excitation energy. In other words, 
generally, F^k/^A^k cx: 1/5 and thus one expects that ef¬ 
fects of decays can be significant only for low-5, quantum 
magnet s.l^ 

Here we offer a general scenario in which this seemingly 
invincible logic fails dramatically. If both decay products 
belong to the flat modes with constant energy 5i, a re¬ 
markably stronger effect must be taking place. Namely, 
in this case the self-energy of the dispersive modes ex¬ 
hibits an essential singularity at the energy 2 ^ 1 , and, for¬ 
mally, the linewidth F^k in 0 is infinite at this energy, 
the effect we refer to as the resonance-like decay. 

In reality, quantum fluctuations of the same origin, i.e. 
coming from the anharmonic cubic terms, a lso ge nerate 
effective further-neighbor J 2 spin couplingswhich 
necessarily warp the flat mode and thus provide natu¬ 
ral means of regularizing this essential singularity. How¬ 
ever, the resultant fluctuation-induced bandwidth of the 
flat mode is now 5-independent, fe^yq oc 0(5^), so that 
the regularized resonance-like broadening in the vicinity 
of 2^1 must scale together with the excitation energy: 
r^k cx^ ^Aik oc 5. This qualitative consideration implies a 
spectacular quantum effect: a very strong damping, elim¬ 
inating spectral weight from the respective energy range 
even in large-5 magnets. Thus, frustration provides nec¬ 
essary and sufficient ingredients for the proliferation of 
intrinsically quantum phenomenon of decays into inher¬ 
ently classical spin systems. 

Altogether, we predict that anomalous broadening and 
a wipe-out of the spectral weight, associated with the 
resonant-like decays, should be common in the spectra 
of the flat-band frustrated systems. In practice, we ar¬ 
gue that the quasiparticle breakdown with characteristic 
termination points and ranges of energies dominated by 
broad continua must be present in 5 = 5/2 kagome-lattice 
Fe-jarosite. 


III. FE-JAROSITE 

In realistic kagome-lattice antiferromagnets, the de¬ 
generacy within the manifold of classical 120° states 
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FIG. 2: (Color online) (a) Neutron-scattering data frorrP^ along the FXYF path (inset). Lines are linear spin-wave theory hts 
of the dispersive modes (dashed) using 0. and the flat mode (solid) with J 2 added to parameters are as shown. Lower 
shaded area highlights the flat band, upper is the set of energies of two flat modes. Arrows imply a decay process into two flat 
modes, (b) Lower curve with the shading is the on-shell Fk from §. Dashed line is the linear spin-wave theory energy of the 
gapless dispersive mode from (a), shaded area shows the half-width boundaries of a lorentzian peak, Sk ± Fk. 


is, most commonly, lifted by the synmetry-breaking 
Dzyaloshinskii-Moriya (DM) terms,I^SES|3Zl yi elding the 
Hamiltonian that cl osely describes Fe-jarosite^^^ and 
other system^^^MEl] 

•H = ^ (JSi • Sj + D • Si X S^) , (4) 

i'ij) 

where summation is over the nearest-neighbor bonds and 
D = (0,0, TDz) on the up/down triangles with the order 
of the site indices in Q shown in Fig. [^b). The out-of¬ 
plane DM interaction lifts the degeneracy and selects the 
q = 0 ground state, i.e. a “ferro”-120° pattern, Fig.j^a). 
A small in-plane DM ternP^ is neglected for simplicity. 

Given the large spin value, S = 5/2, we estimate that 
the ordered moment should be nearly 90% of its classi¬ 
cal value.ESl Similarly, the results of the earlier neutron 
scattering in Fe-jarosit^^ have been inter pre ted as fully 
describable by the linear spin-wave theory,!^ a construc¬ 
tion whose validity we question next. 

Our Fig. I^a) shows the linear spin-wave theory fits 
of the neutron-scattering datgP^ using model Q where 
three distinct excitations branches are easy to identify. 
The DM anisotropy shifts the flat mode from zero en¬ 
ergy to 5ik ~ where du = V^Dz /see Ap¬ 

pendix The flat mode is also not entirely flat. This 
was interpreted^^ as a sign of a phenomenological next- 
nearest-neighborsmer exchange J 2 , ignoring its possible 
quantum origin.^^U gince in the following we do not at¬ 
tempt a fully self-consistent calculation, the same inter¬ 
pretation suffices, with an explicit expression for the dis¬ 
persive flat mode given in Appendix Aside from this 
slight dichotomy with the origin of the flat-mode dis¬ 
persion, the linear spin-wave theory seems to provide a 
spectacular account of the data without the need of any 
quantum effects. 


However, we point out that the spectral weight is con¬ 
spicuously missing from experimental data in the range 
of energies 15 — 19meV in Fig. [^a), i.e. no signal has 
been detected there. While this feature has not been 
emphasized in Ref. [15] and one may argue that the col¬ 
lected experimental data points were simply too sparse, 
the missing band is strongly implied by our discussion, 
as it is exactly in the range of twice the energy of the flat 
mode, 25ik, see Fig. [^a). 

In Fig. I^b) we present the results of the on-shell cal¬ 
culation of Fk for the gapless dispersive mode using § 
with the flat-mode dispersion induced by J 2 for the same 
parameters as in Fig. |^a), see Appendix [A| for details. 
As we discuss later, the dynamical structure factor allows 
to view modes selectively in different parts of the k-space 
and in different polarizations.^SlThe results for the damp¬ 
ing are combined with the energy s\^ of the mode with 
the shaded area representing half-width boundaries of a 
lorentzian peak, ^kiTk- We have also verified^ that the 
effect of renormalization on the real part of the spectrum 
is minor, in agreement with approximation in 

Our Fig. [^b) demonstrates that the spin-wave exci¬ 
tation is well-defined until a sharp threshold at about 
25^^. Above that energy, the broadening reaches about 
one third of the bandwidth signifying an overdamped 
spectrum, consistent with the missing spectral weight 
in the experimental data. The sharp transition implies 
a threshold singularity and other spectral features that 
are characteristic to the quasiparticle breakdown phe- 
nomen on in quantum Bose liquids and S = 1/2 spin- 
liquids.l^^Sl] xhere is a partial reconstruction of the spec¬ 
trum at the energies above 2^^^ where decays are no 
more resonant-like as indicated in the figure, i.e., occur¬ 
ring due to other, non-resonant channels, but still pro¬ 
viding a sizable broadening to the spectrum. 
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FIG. 3: (Color online) (a) A 3D plot of the magnon dispersion for Fe-jarosite within the linear spin-wave theory with planes 
of the cuts in (b) and (e). (b)-(f) Intensity maps of Aq(ct;) in units of {2SJ)~^ vs q throughout the Brillouin zone for a set of 
energies. Intensity scale is as described in text. Dashed black lines are peak positions from the linear spin-wave theory. 


The non-resonant decays result in a typical broad¬ 
ening F ~ 0.25J, i n acc ord with similar results for 
the triangular-latticj^^^ and other frustrated spin 
systems.^ By contrast, the broadening in the resonant- 
decay region in Fig.j^b) reaches F?^ 1.7J, an effect larger 
by a factor exceeding 2S for the considered 5' = 5/2 model 
of Fe-jarosite. This is in a remarkable agreement with our 
qualitative discussion on the scaling of the resonance-like 
decay rate with S', provided after Eq. (§ above. 

We note that the broadening on the top of the band in 
the non-resonant region translates to less than ImeV, be¬ 
low the experimental resolution of Ref.[T5]in which all the 
data were described as resolution-limited. The current 
resolution of the neutron-scattering experiments is easily 
an order of magnitude higher. We also point out that our 
consideration is aimed at the strong qualitative features 
of the spectrum of a representative flat-band frustrated 
spin system, not on the minor quantitative details. As 
such, small discrepancies with some of the data may oc¬ 
cur due to, e.g., neglect of the in-plane DM terms, but 
should be considered as secondary. 


Dynamical structure factor 

To demonstrate the effect of decays, we performed a 
calculation of the magnon spectral functions, Ajyq(u;), 
quantities directly related to the spin-spin dynamical cor¬ 


relation function via 

>S““(q,w) (X f (5) 

U 

Here, the kinematic formfactors allow to “filter out” 
spectral contributions of some of the modes to the in¬ 
plane and the out-of plane components of tS(q, cj) in the 
portions of the q-space while highlighting the other ones: 
a phenomenon akin to the extinction of the Bragg peaks 
in the non-Bravias lattices.l^ Using this feature, we con¬ 
centrate only on one of the dispersive modes. 

A dramatic view on the drastic transformations of the 
spectrum can be observed in constant-energy cuts of the 
dynamical structure factor in the range of energies af¬ 
fected by the resonance-like decays. In Fig.|^ we present 
intensity maps of such constant-energy cuts for Ajyq(u;), a 
close proxy of the dynamical structure factor 5(q, cj), for 
the dispersive magnon mode for the energies ranging from 
Il.TmeV to 20meV. The upper cut-off of the spectral 
function is chosen to correspond to the maximal height 
of the peaks in the non-resonant decay region in Fig.|^b) 
and translates into the broadening Fk~0.73meV for the 
Fe-jarosite values of S and J, which should be resolvable 
by the modern neutron-scattering measurements. 

The first of the cuts is below the threshold energy 
2^5^^ and shows a very close accord of the sharp-intensity 
peaks in with the expectations from the linear, 

non-interacting spin-wave theory, shown by the dashed 
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lines. The three subsequent cuts, Fig. [^c)-(e), are from 
within the resonant-decay band, 2^5^^ <uj< where 

one can observe strong deviation from such expectations, 
massive redistribution of the spectral weight into differ¬ 
ent regions of the q-space, and a multitude of intriguing 
“shadow” features, reflecting van Hove singularities in 
the two-particle density of states of the decay products.^!! 
The last cut. Fig. if), is, nominally, above the top of the 
magnon band and should be expected to show zero inten¬ 
sity everywhere. Instead, it is also affected by the spec¬ 
tral weight redistribution and retains some of the features 
of the other cuts. Altogether, Fig. [^b) and Fig. [goffer 
a comprehensive theoretical insight into the non-trivial 
features of the dynamical structure factor of a flat-band 
kagome-lattice antiferromagnet, which originate from the 
decays of magnetic excitations facilitated by the nonlin¬ 
ear couplings. 


IV. SUMMARY 

To summarize, we have outlined a general scenario for 
drastic transformations in the spectra of frustrated mag¬ 
nets that feature flat modes and have substantiated it 
by a consideration of the spin-spin structure factor of 
the large-5' kagome-lattice system Fe-jarosite. Our study 
calls for further studies in these systems. 

We would also like to comment that recently, the broad 
features in the spectra of magnetic systems have become 
a direct sign of fractionalized excitations of prospective 
spin-liquid phases.l^^®^ In this work, we have provided 
a case study of an excitation spectrum of a strongly 
frustrated but almost classical and well-ordered kagome- 
lattice antiferromagnet, for which we have demonstrated 
extremely strong broadening and even a complete and 
spectacular wipe-out of a part of its spectrum. Here, the 
broad features are due to flat or weakly dispersive modes, 
a hallmark feature of a variety of frustrated spin systems, 
and due to a non-collinearity of spins in the ground state, 
again an outcome of competing interactions. Thus, this is 
also a cautionary tale, because the same reasons that may 
lead to the spin-liquid behavior may also favor strong 
coupling and decays among quasiparticles. 
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Appendix A: Technical details 
1 . Spin-wave theory 

Eollowing the approach outlined in Refs. 181121 one can 
diagonalize the harmonic part of the Hamiltonian in Q 
to obtain the spin-wave energies 

5ik = 2JS\/3dM (1 + ^m) /2 , (Al) 

for the “flat mode,” and 

^2(3)k = 2 J5\/l -h dM (^2) 

X + dM — 7k — dM (1 ± \/l + 87 k )74 , 

for the dispersive modes, dM = V^Dz/J here. 

Corrections due to effective J 2 interactions can be 
taken into account perturbatively to yielcP^ the disper¬ 
sion of the “flat mode” 

£ik = 275^(3(1+ dM )/2 + j 2 (i-A(^^/ 2 )) (A 3 ) 
{^M + 32 (1 + ) , 

where ^2 = J 2 IJ and 

= (/2(k) - /i(k)) / (1 - 7k) , (A 4 ) 

with /i(k) = 4ci + C2C2 + C3C3 , 

/2(k) = C1C2C3 + C2C1C3 + C3C1C2 , 

with the shorthand notations = cos{qn), c[=cos{qs-\- 
^ 2 ), c '2 = cos(g 3 -gi), ^ = cos(gi+g 2 ), where gn = k-^n/2, 
and Sn are the primitive vectors of the kagome lattice. 

The diagonalization of the harmonic part of Q implies 
a two-step procedurd^^ with the unitary transformation 
of the original Holstein-Primakoff bosons 

^ak ~ ^ ^ 'a)jy^Q;(k) djyk 5 

V 

followed by the usual Bogolyubov transformation for 
each of the individual species of d-boson, see Ref. [12] 
for details and for the explicit form of the eigenvectors 
= (w^,i(k),M;^,2(k),M;,.,3(k)). 

2 . Cubic terms 

Due to noncollinear 120° spin struct ure, c ubic anhar- 
monic coupling of the spin waves occurs.l^^^^It originates 
from the 5f 5| terms in Q, written in the local reference 
frame.^ In the bosonic representation they yield 

% = J (1 + dMsin 0 ij{ala]aj +h.c.) , (A 6 ) 
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where Oij = ±120° is the angle between two neighboring 
spins. 

Assuming the spins in the q = 0 state, Fig. [^a), and 
using the unitary and Bogolyubov transformations men¬ 
tioned above gives the “source”, and the “decay”, 

6^6^6, terms, see Ref. [12] where the effects of the former 
were discussed. The decay part of the Hamiltonian is 

± = E ^qk?p±±^»7P+h.c., (A7) 

' ^ ^ k+q=p 


with the vertex 


qk;p 



qk;p ■ 


(A8) 


for the resonance-like decay phenomena, one can approx¬ 
imate the self-energy by its on-shell imaginary part, i.e. 


~ RmE^k(^/xk) — '^F^k . 


(A12) 


which is given by Clearly, the dispersion of the flat 
mode in (A3) is crucial for the decays into two of them. 


as otherwis e this channel would produce essential sin¬ 
gularity in (All) and in F^k- With that, evaluation of 
the spectral function A^\^{uj) = — (l/7r)ImG^k(<^) can be 
performed numerically. 

The diagonal components of the dynamical struc¬ 
ture factor, or the spin-spin dynamical correlation func¬ 
tion, which contribute directly to the inelastic neutron¬ 
scattering cross section, are given by 


which is explicitly oc ^/S. The symmetrized dimension¬ 
less vertex ^qk^p is given by 



dt 

2tt 




(A13) 


^qklp “ “I” '^vq){u^^Urip + Vfik^rip) 

+ ^kpq(“Atk + V^k)( ± ^zyp^?7q ) (A9) 

± -^pqk(^’7P "Ct^p) ('^z^q'^^ik ± '^zyq'^/xk): 

where Uy\^ and Vy\^ are the Bogolyubov parameters and 
the amplitudes Tqkp given by 


al3 


where is the Levi-Civita antisymmetric tensor, and 
shorthand notations are = qp/3a and = Pp — Pa^ 
here pa are the atom’s positions within the unit cell. 


3. Self-energy, spectral function, and structure 
factor 


Using the standard diagrammatic rules for ([^, we ob¬ 
tain the second-order decay self-energy 






up 

I q,k—q;k I 


^ ^T] ^ ^77k-q + 


(All) 


which contributes to the l/5'-correction to the magnon 
energy. The magnon Green’s function for the branch p is 
given by (§• Since only the decay terms are responsible 


where ao refers to the laboratory frame {xo,^ 0 :^o}- 
Given the co-planar spin configuration, it is convenient 
to separate the in-plane and out-of-plane components 
of 5^°^(q, cj). Assuming equal contribution of all three 
ao components to the cross section, using the spin-wave 
mapping of spins on bosons with the two-step transfor¬ 
mation described above, after some algebra, one can ob¬ 
tain the leading contributions to the structure factor as 
directly related to the spectral function 

5in(out) A,q(a;), (A14) 


where are the kinematic formfactors. It is im- 

portant to note that the kinematic formfactors are mod¬ 
ulated in the q-space and are suppressed in one of the 
Brillouin zones while are maximal in the others. ESI This 
effect is characteristic to the non-Bravias lattices and is 
similar to the effect of extinction of some of the Bragg 
peaks in them. Because of that, one may be able to high¬ 
light spectral contribution of one of the magnon branches 
while “filtering out” the others by selecting a particular 
component of the structure factor in a particular Bril¬ 
louin zone. Our analysis demonstrates that the out-of 
plane component of tS(q, cj) should be totally dominated 
by only one of the dispersive modes (gapless) in one of 
the three distinct Brillouin zones. This feature can be 
useful for the future neutron-scattering experiments. 
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